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Abstract. Predictions of analytical models for hadronisation, namely the dispersive model, the shape 
function and the single dressed gluon approximation, are compared with moments of hadronic event shape 
distributions measured in e^e^ annihilation at centre-of-mass energies between 14 and 209 GeV. In contrast 
to Monte Carlo models for hadronisation, analytical models require to adjust only two universal parameters, 
the strong coupling and a second quantity parametrising nonperturbative corrections. The extracted values 
of Qs are consistent with the world average and competitive with previous measurements. The variance of 
event shape distributions is compared with predictions given by some of these models. Limitations of the 
models, probably due to unknown higher order corrections, are demonstrated and discussed. 

PACS. 12.38.Lg Other nonperturbative calculations - 12.38.Qk Experimental tests 



1 Introduction 

In previous studies, moments of event shape distributions 
have been compared to perturbative predictions of quan- 
tum chromodynamics (QCD, [1-4]) in next-to-leading or- 
der, simulating the hadronisation process by Monte Carlo 
models [5,6]. Alternatively there exist models describing 
hadronisation analytically. This paper aims to study these 
models qualitatively and quantitatively by measurements 
of the strong coupling and the model parameters. 

The dispersive model [7] is based on the assumption of 
a nonpertubatively continued strong coupling. The shape 
function [8] additionally describes a modification of the 
shape of the perturbatively calculated distribution. The 
single dressed gluon approximation [9, 10] estimates the 
perturbative part more completely with reduced pertur- 
bative uncertainties of the prediction. The models test the 
predicted energy evolution of the strong coupling. Their 
parameters can be determined — one of them the value of 
as at some reference energy — and the assumption of uni- 
versality of these parameters can be probed. To date, pri- 
marily the distributions themselves and the mean values 
(first moments) have been studied [11-14]. 

This analysis uses data measured by JADE [5] in the 
years 1979. ..1986 at six centre-of-mass (cms) energies in 
the energy range of Q = 14. ..44 GeV, and data measured 
by OPAL at 12 energy points over the whole LEP en- 
ergy range of 91... 209 GeV and combined into 4 energy 
ranges [6]. The large energy range covered by the mea- 
surements allows to test the employed assumptions selec- 
tively. 



The outline of the paper is as follows. In Sect. [U we 
present the observables used in the analysis. In Sect. [3] 
we describe the perturbative QCD predictions and intro- 
duce the analytical models which describe the hadronisa- 
tion process. In Sect. [4] we discuss predictions for event 
shape moments at hadron level and compare them with 
the measurements. The predictions for the event shape 
distribution variance following from the dispersive model 
and the shape function are tested as well. In Sect. [5] we 
summarize and give our conclusions. 



2 Event shape moments 

Event shape variables are a convenient way to characterise 
properties of hadronic final states. They are calculated 
from particle momenta and energies. For definitions of the 
variables we refer to [6]. 

In a hadronic event in e+e" annihilation the virtual 
vector boson 2^/7* generated by annihilation of electron 
and positron decays into a quark pair qq. The quarks may 
radiate gluons which radiate further gluons or decay into 
another quark pair. The final state of this parton shower 
is called parton level. By the process of hadronisation par- 
tons are transferred into hadrons. The predictions used in 
this work describe the hadron level. The variables mea- 
sured in the experiments have to be corrected for the 
effects of limited detector acceptance and resolution to 
probe the hadron level. 

The event shapes considered here are thrust T, C- 
parameter C, heavy jet mass Afu, jet broadening vari- 
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ablcs i?T and -Bwj a-nd the transition value j/^ between 
2- and 3-jet final states defined using the Durham jet 
reconstruction scheme. The as determination in [5] is 
based on moments of these variables, in [6] on distribu- 
tions and moments of these variables. Their theoretical 
description by perturbation theory is currently the most 
advanced [15-17]. 

A generic event shape variable is denoted by the sym- 
bol y. Regions dominated by multiple jets give large val- 
ues of y, while two narrow jets give ?/ ~ 0. Thrust T is 
an exception to this rule. By using y ~ \ — T instead, 
the condition is fulfilled for all event shapes. i3vv, 2/23 and 
Mh are sensitive to only one suitably chosen hemisphere 
of the event [one-hemisphere variables) , 1 — T, C and i?T 
are sensitive to the whole hadronic event [two-hemisphere 
variables) . 

The nth, n = 1,2,... moment of the distribution of an 
event shape variable y is defined by 



1 dcr 



(1) 



where j/max is the kincmatically allowed upper limit of the 
variable and crtot. denotes the total hadronic cross section. 

Predictions have been made available for the moments 
of event shapes. Their evolution with cms energy allows 
direct tests of the predicted energy evolution of the strong 
coupling as . Furthermore it enables the determination of a 
single value of as at a definite energy scale — for example 
as(M2o) at the rest energy of the Z° boson. The theo- 
retical calculations always involve a full integration over 
phase space, which implies that comparison with data al- 
ways probes all of the available phase space. This is in 
contrast to QCD predictions for distributions; these are 
commonly only compared with data — e.g. in order to mea- 
sure as — in restricted regions, where the theory is able to 
describe the data well, see e.g. [18]. Comparisons of QCD 
predictions for moments of event shape distributions with 
data are thus complementary to tests of the theory using 
distributions. 



3 Theory 

The QCD prediction for (y") in next-to-leading-ordeir 
(NLO) of the strong coupling as = as/(27r) has the forrrO 



(y") = An 



as + [Bn - 2An)al . 



(2) 



^ Very recently, perturbative NNLO predictions for event 
shape moments became available [19]. The nonperturbative 
models discussed in the following sections need to be adapted 
to these new NNLO predictions before a meaningful study 
based on NNLO can be done. Therefore, here we restrict our- 
selves to consistently using NLO together with the currently 
available models. 

^ The a| coefficient is written as B„ — 2An because the QCD 
calculations are normalized to the Born cross section ctq, while 
the data are normalized to the total hadronic cross section 
cTtot = o"o(l + 2as) in LO. 



The coefficients An and Bn arc obtained by numerical 
integration of the QCD matrix elements using the program 
EVENT2 [20]. These predictions were also used in [5,6]. 

The coupling as, and therefore the QCD prediction 
depends on the renormalisation scale ^jj, see e.g. [21]. 
The prediction is changed by this dependence as shown 
in [5]. For clearer notation the renormalisation scale fac- 
tor is defined as a;^ = hr/Q] a truncated fixed order QCD 
calculation such as ([2]) will then depend on x^. 

Infrared and coUinear stability are essential for a per- 
turbative description of the parton level. These proper- 
ties, however, do not suffice for a perturbative description 
of the hadron level as hadronisation takes place at low 
energy scales, where the perturbative description breaks 
down. The evolution of partons to hadrons can be approx- 
imated by analytical calculations. These analytical calcu- 
lations are generally motivated from the transition from 
perturbative to non-perturbative regime. A perturbation 
series in quantum field theory is a divergent series, see 
e.g. [22]. To obtain finite results, regularisations have to 
be applied. From suitable prescriptions non-perturbative 
terms are found, which typically scale with inverse powers 
of the cms energy [7-9]. 



3.1 Dispersive model 

This model [7,23,24] is based on the assumption of a non- 
perturbatively defined strong coupling as(Q^) which re- 
mains finite at and below the Landau pole vIqcd- The Lan- 
dau pole is the scale where the usual perturbative coupling 
diverges. The matching scale Q = , marking the border 
between perturbative and non-perturbative region, is not 
uniquely defined. Usually it is taken as yU/ ~ 2 GeV. 

As the non-perturbative coupling cannot be calcu- 
lated, it is parametrised universally in a simple way by 
the zeroth moment of the extended coupling over the non- 
perturbative region. 



1 f^' 
ao[ni) = — / dQas 
l^i Jo 



(3) 



In first approximation, non-perturbative corrections 
generate a simple shift of the perturbative differential dis- 
tribution dffj^^/dy of the event shape variables 1 — T, C, 
St, anq^A/H when relating parton to hadron level. 



dy 



dy 



[y-Oy- V) 



(4) 



This prediction is valid only if the value of the event shape 
variable y is not too large (y <IC 1), and the cms energy Q 
not too small [Q Aqcd/u)- Only the numerical factor 
Oy depends on the event shape variable y, see Table [T] 
However V depends on the hard scale Q but is universal 



The theoretical calculations are based on the variable Mh 
because of the problems with the NLO description of Mh as 
discussed in [5,6,25]. 
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for the event shape variables 1 — T, C, Bt , -Bw and AI^, 
and has the form [23,24] 



P 



ACf 



+2/3oa^(Ai^) ln^ + 1 



(5) 
/if 

Q 



In the MS renormahsation scheme the constant K has the 
value 

K^cJ'l-^-\-'-n, (6) 

Vi8 6 y 9 ^ ' 

with rif = 5 at the studied energies, and the beta func- 
tion coefficient is /3o = 23/(127r). The color factors have 
the values Cp = 4/3, Ca = 3 [26], the Milan factor M is 
known in two loops, A4 = 1.49 ± 20% (for flavour number 
rif = 3 at the relevant low scales). The cited uncertainty 
results [24] from the estimation of the next order contri- 



bution, tVI, 



(l + 0(asA)). 



Table 1. Coefficients ay of power correction oc 1/Q of event 



shape variables in 


the dispersive 


model 


[27,28] 




event shape 
variable y 


1-T C 


Bt 


Bw 


V2Z Ml 


ay 


2 37r 


1 


1/2 


1 



Applying the dispersive model calculation ^ of the 
normalized event shape distribution in definition ([1]) of the 
moment of order n and naively neglecting the integration 
over the unphysical range of negative variable values, gives 



1 1 

ay 





dy {y + ayVr 



dy 



-iy)- 



The predictions for the moments on hadron level become: 

{y') = {y')..o+ayr (7) 

(y') = (?/')nlo + 2(y^)„,o ■ayV + (ay rf (8) 
+ {ayVf (9) 

(y') = (y')NLo + 4(y'),,o ■ayV + 6{y')^^o ■ (dy rf 
+4(y')..o-K^)' + (s^)' (10) 

{y') = {y\.o + 5(y')«.o ■ «y ^ + io(y')«.o • K 
+io{y^)^^a ■ (% + Hy')..o ■ K 
+K^)'- (11) 

Previous studies [29] indicate that the parameters 
as{Mzo) and Q;o(/i/) when fitted to and _Bw distribu- 
tions via ^ are not compatible with values derived from 
1 — T and C. Therefore improved predictions for these dis- 
tributions were given. They also describe a compression of 



the distribution peak, i.e. a narrowing in the two jet re- 
gion [24,30]. The non-perturbative factor P in the case of 
{B^) and {B^) is replaced by [31] 



B\) 



2.JCFa,{l + Ka,/{2^)) 

3 27r/3o 
+ 4- ' 



3Cf 



(12) 



rsp. 



3 _ 7r/3o 
4 



3Cf 
«s(c-3/2 



'70 



(13) 



)^) and a 



with a rescaled coupling as{Q^) 
constant ryg — —0.6137. 

No power correction coefficient has been calculated in 
the dispersive model for the variable y^^. The purely per- 
turbative prediction describes the first moments of y^g 
well [14,31] — therefore we also compare it with the higher 
moments. 

The dispersive model gives predictions for several ob- 
servables and contains only universal free parameters 
as(Mzo) and ao(M/)- 



3.2 Shape function 

Korchemsky and Tafat [8] describe properties of the event 
shape variables 1 — T, C and not included in NLO 
perturbation theory by a so called shape function, which 
does not depend on the variable nor the cms energy. This 
is more general than the dispersive model, as not only a 
shift of the perturbative prediction is predicted but also a 
compression of the distribution peak. 

The prediction is deduced from studying the two jet 
region (i.e. y <C 1) in the distribution of the event shape 
variable y. The prediction for the differential distribution 
is 



1 da{Q) 
a- dy 



'^'y da 
def.ie/-^ 



{y - e/Q) , 



(14) 



with a non-pcrturbativc function fy{e), dependent on one 
scale parameter e. This function is derived from the shape 
function f{eL, £r) [8], which depends on two scale param- 
eters £l: £r for the two hemispheres of the event. By the 
compression of the distribution the validity of the pre- 
diction is extended compared to the dispersive model to 
y - ^qcd/Q- 

The shape function can be parametrised [8] by its "first 
moment" , 

Al = J d£R J deLi£B. + £L)f{£R,£L) = {£b. + £l) , (15) 
its "second moment" , 

\2 = {i£R + £L?), (16) 
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and a Q dependent function 5X2{Q). Predictions for the 
moments of event shape variables can be derived, and Ai 
and A2 can be fitted to the data. 

From prediction p4|) for the distribution of the vari- 
ables 1 — r, C and M^, predictions for the mean values 
follow [8] by integration, 



((i-r)i) = ((i-T)i)_^ + 



Al 
Q 



1 - 5.73 



27r 



2Q 



Analogously for the second moments one finds 



{{l-Tf) 



((1-T)2) +2^(1 -T) + 



(02 



(17) 
(18) 
(19) 

], 

(20) 



+ 2 Q 

97r2 A2 



2(C) -4.30 4^ 



4 Q2 



1 - 11.46 

4^ (^4) 



«s(g' 



27r 



(21) 



X2 + SX2iQ) 
4Q2 



(22) 



In this model the more strongly suppressed power correc- 
tions have an independent coefficient. The coefficient Ai is 
interpreted as the first "moment" of the shape function, 
A2 the second, and <5A2 as a contribution accounting for 
the one-hemisphere character of M^. Therefore these are 
universal scales. 



3.3 Variance of event shape distributions 

In [32] the variance of event shape distributions on hadron 
level was found to be described in the dispersive model 
perturbatively without significant power corrections. This 
would open the possibility of an accurate ag {Mzo ) deter- 
mination. 

A theory predicting first and second moments of an 
event shape variable y also gives a prediction for the vari- 
ance Var(y) = (y^) — (y)^, and so a simple prediction for 
the variance of event shape variables on hadron level is 
deduced [32] by employing the predictions ([T]) and 



Var(y) 



(y^)NLO 



{y) 



2 

NLO 



(23) 



In the dispersive model one obtains a purely perturba- 
tive expression for the variance, up to strongly suppressed 
corrections 0{as/Q'^)- 

Predictions for the event shape variance can also be 
derived from the shape function. The first and second mo- 
ment predictions (fT7|) and (|20p give the identical predic- 
tion in case of thrust, 



Var(l-r) = ((l-r)2)^ 



{i-ry 



Analogously this follows from ((T9|) and ([22]) for N&t{M'^. 
Var(C) follows from (HH]), jSlI) and the perturbative LO 
coefficient [20]; up to ©(l/Q^) 

Var(C) = (C^)^^^ - {C)l^^ - 3.23^a.(g2) , (24) 

showing an additional term (x a^/Q with coefficient Ai as 
in Subsect. 13.21 



3.4 Single dressed gluon approximation 

Gardi et al. [9, 10] assume the existence of a reordering 
of the perturbative series, the so called skeleton expan- 
sion [33] — its existence is proven only for abelian field the- 
ory like QED. The first contribution to this expansion is a 
single dressed gluon (SDG) which resums running coupling 
effects of any order in as- These are renormalons in the 
dominant contributions oc /3q . The single dressed gluon 
graphs can be calculated completely. In this way the per- 
turbative prediction ((1 — T)")pt. for a moment of thrust 
can be approximated. This approximated series already di- 
verges and various regularisations differ by definite powers 
of the cms energy Q. These predictions follow: 



((1 




= (1- 


-T) +^ 


((1 




= ((1 




((1 


-Tf) 


= ((1 




((1 


-Tf) 


= ((1 





V2 


















Ka 







In general, the non-perturbative correction is predicted 
as a sum of two terms with different powers in Q and 
coefficients i^n, k„. 

The thrust is used in the massless limit. 



T 



(25) 



Here the denominator is changed with respect to the stan- 
dard definition [6] by X^i X^i-^*' which does not 
change the thrust value as long as massless partons are 
considered. Under this change of definition the thrust 
value calculated with a "massive" virtual gluon is cor- 
rect, as long as all (massless) partons generated in the 
hadronisation end up in the same event hemisphere w.r.t. 
to the thrust axis. This so called "inclusive" calculation is 
valid only if fragmentation of the gluon is approximately 
coUinear. The perturbative SDG predictions for the mo- 
ments of thrust are given as perturbation series in the 
coupling a [10], which relates to as{f^) in the MS scheme 
by a simple shift of the Landau PolcQ, 



1 



(26) 



* In particular the beta function coefficients /3i have the same 
values as in the MS scheme. 
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The pcrturbative cocfRcicnts arc given in so called log 
moments di of the characteristic thrust function. The first 
six log moments for ((1 — r)^)...((l — T)^) were calculated 
by numerical integration [10]. The pcrturbative prediction 
of C(ag) in the (obviously incomplete) SDG approxima- 
tion has the form [34] 



= do • a + /3o dxTx ■ d^ 

+ ({-Ido + d2)l3l + pi di^ TT^ ■ d^ 

+ (^{-7T^ di + d3)Pl + ^[3iPod2+ (32 d^ TT^ • 

+ (^{d4 + ^doTT*-2n^ d2)f3^ + (y /3i d3 - I3i tt^ di)l3l 

+3 132 d2 /3o +^ Pi d2 + /?3 di) TT^ . 



/ 10 77 q 



-(6 da - di)f52 Pi + (y Pi d, + -P, d2)Po 



+/34di + -/3i/32d2 )^'-a'; 



(27) 



where the n-dependent di are taken from [10]. Higher or- 
der contributions from the running of the coupling are 
accounted for up to sixth order. Reference [10] gives the 
prediction explicitly in 0{d^) without the term Pidiw'^ -d^ . 
This term with marked effect on the a^-coefficicnt re- 
sults from considering the usual QCD beta function [34]. 
The symbols P0...P3 denote the usual beta function coef- 
ficients [35,36], P4 is unknown. We set P4 to 0. 

The SDG approximation (P7)) is complete in leading 
order 0{d) by construction [10]. In higher orders this ap- 
proximation only gives the terms oc Po- Therefore we use 
the numerically calculated [20] NLO coefficients. For the 
third and higher orders we employ the SDG approxima- 
tion. 



(28) 



Dn -d^ + En-d^ +F„ ■ d^. 



As the pcrturbative expansion is an asymptotic (diver- 
gent) series [22], the terms are expected to become smaller 
up to a certain order, and then become larger again. It is 
of interest 



— at which expansion order this happens, 

— how the measured coupling depends on the maximum 
expansion order, 

— how the leading power correction i-'„/(5'" depends on 
the maximum expansion order, 

— how the minimum pcrturbative term relates to the 
power correction. 

The best approximation of the theory by an asymptotic 
series is expected when truncating it near the minimal 
term — including additional terms does not necessarily re- 
sult in a better approximation. Therefore we study the 
analysis as a function of the truncation order. 
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4 Tests of non-perturbative models 

Our tests use event shape moments with statistical and 
experimental uncertainties measured by JADE and pub- 
lished in [5] and analogous OPAL data published in [6], 
covering in total the energy range of 14 to 209 GeV. In the 
JADE data e"''e~ — > bb events have been subtracted on 
a statistical basis. Table [2] gives an overview of the data 
used. 



Table 2. Year of data taking, energy range, integrated lumi- 
nosity, average cms energy and the numbers of selected data 
events for each JADE [5] or OPAL [6] data sample. The hori- 
zontal lines in the OPAL ranges separate the data into the four 
energy ranges used for fit and presentation purposes 



year 


range of 
Q in GeV 


mean Q 
in GcV 


luminosity 
(pb-1) 


selected 
events 


1981 


13.0 ., 


,. 15.0 


14.0 


1.46 


1783 


1981 


21.0 ., 


.. 23.0 


22.0 


2.41 


1403 


1981, 1982 


33.8 ., 


,. 36.0 


34.6 


61.7 


14313 


1986 


34.0 ., 


,. 36.0 


35.0 


92.3 


20876 


1985 


37.3 ., 


,. 39.3 


38.3 


8.28 


1585 


1984, 1985 


43.4 ., 


,. 46.4 


43.8 


28.8 


4376 " 


1996, 2000 


91.0 ., 


.. 91.5 


91.3 


14.7 


395695 ' 


1995, 1997 


129.9 ., 


.. 136.3 


133.1 


11.26 


630 ' 


1996 


161.2 ., 


.. 161.6 


161.3 


10.06 


281 " 


1996 


170.2 ., 


,. 172.5 


172.1 


10.38 


218 ' 


1997 


180.8 ., 


,. 184.2 


182.7 


57.72 


1077 ' 


1998 


188.3 ., 


,. 189.1 


188.6 


185.2 


3086 ' 


1999 


191.4 ., 


.. 192.1 


191.6 


29.53 


514 1 


1999 


195.4 ., 


.. 196.1 


195.5 


76.67 


1137 ' 


1999, 2000 


199.1 ., 


,. 200.2 


199.5 


79.27 


1090 ' 


1999, 2000 


201.3 ., 


,. 202.1 


201.6 


37.75 


519 1 


2000 


202.5 ., 


,. 205.5 


204.9 


82.01 


1130 1 


2000 


205.5 ., 


,. 208.9 


206.6 


138.8 


1717 1 



JADE and OPAL are similar in construction and 
many parameters [37]. Consistent measurements can be 
expected from the simultaneous use. The analysis pro- 
cedures for both data sets were constructed to be simi- 
lar. Systematic variations of the JADE and OPAL analy- 
ses concern detector event reconstruction, selection cuts, 
Monte Carlo generators and background. 



4.1 Dispersive model 

Our test uses the first five moments of 1 — T, Bt , i?w : 
C and j/^, and the second and fourth moment of A/h- We 
compare the prediction with the hadron level data, varying 
two parameters asiM^o) and ao{fJ.i)- Figs. [H [2] show the 
comparison of data and predictions ([7l. fTT|) in the case of 
1 — T and _Bw The x^/d.o.f. values — calculated from the 
statistical errors only — vary between 2 and 10, but should 
be regarded as being indicative only. Experimental differ- 
ences between JADE and OPAL contribute significantly 
to the x^/d.o.f . values. At low cms energies the predictions 



drop off again — this unphysical behaviour (not always vis- 
ible in the figures) substantially contributes to the high 
X^/d.o.f. valuesO The moments of the one-hemisphere 
variables Bw and in particular Mh are described better 
than in the respective comparison with hadronisation cor- 
rection by Monte Carlo models [5] . 

To estimate the experimental systematic uncertainties, 
the fits are repeated based on the minimum overlap as- 
sumption \5\ for combining the OPAL and JADE system- 
atic errorsH 

We study uncertainties of theory parameters by the 
following variations 

— The renormalisation scale uncertainty (i.e. the effect 
of higher perturbative orders) is estimated by setting 
Xfj, = 0.5 and = 2.0. Our default is always = 1. 
The hadronisation uncertainty is estimated by setting 

- M = 1.49 ±20%, 

— — 1 GeV and /i/ = 3 GeV. The resulting devi- 
ations are included into the as{Mzo) uncertainty 
only, as ao(M/) depends on /x/ directly by its defi- 
nition ([3]). 

Tables [3] and |4] and Fig. [3] contain the as (M^o ) and 
ao{fij) results from the standard measurement of the mo- 
ments of the six event shape variable^ and from the 
■a;^. III and A4 variations. The values of as(-Mzo) and 
measured from ((l-T)!}, (C^), (B),), ((yi^a)'), 
(Afy) are in agreement within errors with previous analy- 
ses, sec [31,39] and references therein. The as{Mza) value 
from (B^) is lower and the ao(/i/) value higher than in 
previous analyses because of the incomplete description of 
B^) in NLO discussed below. This incompleteness affects 
the predictions at low centre-of-mass energies, leading to 
high ao{pi) values and low as values. Low energy data 
have not yet been analysed with a statistics comparable 
to ours. 

The as{Mzo) results arc similar to those in the analysis 
employing Monte Carlo models for hadronisiation correc- 
tion [5,6]. In particular the asiM^a) values steeply rising 
with moment order for the variables 1 — T, C and Bt 
show the incomplete description of the multi jet region in 
NLO as already seen and discussed in [5,6]. This is also 
seen in the large renormalisation scale uncertainties. The 
aa{pi) and as(Afzo) values from the moments (B^) are 
not universal but decrease with moment order n. We can 
validate only the cxplicitely calculated (i?^^) prediction. 

Based on the statistical errors, the values of ao{i^i) 
from the (B^), (A/^) and (Af^) fits are significantly 



'' The case of j/^ where we compared with the purely per- 
turbative NLO prediction is different. The energy evolution of 
the data appears consistent with the prediction but the mea- 
sured higher moments at 14 GeV are too low. A suitable power 
correction would scale with a higher power of cms energy, and 
indeed power terms In Q/Q^ and are expected [7] and 

describe distributions better [38]. 

^ The x^/d.o.f. values from these fits are substantially lower 
than those from the fits employing the statistical errors. 

^ We list only results from moments that do not show prob- 
lems in the multi jet region dicussed below. 
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Fig. 1. Fits of the dispersive prediction to JADE and OPAL 
measurements of 1 — T moments. The solid line shows the pre- 
diction with fitted values of asiM^o) and ao{fJ^i), the dashed 
line shows the pure NLO contribution. The inner error bars 
show the statistical uncertainties used in the fit, and the outer 
error bars show the combined statistical and experimental sys- 
tematic errors 
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Fig. 2. Fits of the dispersive prediction to JADE and OPAL 
measurements of Bw moments. The solid line shows the pre- 
diction with fitted values of asiM^^o) and ao{fJ,i), the dashed 
line shows the pure NLO contribution. The inner error bars 
show the statistical uncertainties used in the fit, and the outer 
error bars show the combined statistical and experimental sys- 
tematic errors 



higher than those from ((1 -T)i), (C^) and (B].). The 
NLO description of i3w and Mh in the two jet region 
is rather incomplete compared to 1 — T, C and Bt — 
especially at low energy where the value of the coupling 
is high. At low Q the NLO predictions of the Bw, 2/23 
and Mn distributions are (unphysically) negative in a 
large range of the two jet region [25,38]. This incom- 
pleteness of the NLO prediction for the moments — in the 
case of {B^) the Og coefficient is even negative — is com- 
pensated by a larger power correction. The perturbative 
description then contributes less as is seen by the low 
as{Mzo) values. Fits to JADE data of (B^) alone give 
as{Mzo) = 0.1029 ± 0.0016, ao(M/) = 0.728 ± 0.017; to 
OPAL alone they return as{Mzo) = 0.1242 ± 0.0025, 
ao{fj.i) = 0.234 ±0.103 (statistical errors). So at low en- 
ergies where the coupling is large, the compensation is 
stronger, and the power correction is not universal (JADE 
and OPAL agree better on every other moment). 

For averaging as{M-z,o) and ao{ni) over many variables 
and moment orders, we exclude moments that suffer from 
the problems discussed in the preceding two paragraphs: 
The deficiencies in the description of the multi jet region 
for the two-hemisphere variables have already been seen 
in [5,6]. To select observables with an apparently converg- 
ing perturbative prediction, we consider as in [5, 6] only 
those results for which the NLO term in equation ^ is 
less than half the corresponding LO term (i.e. \K\ < 25 



or \Kas/2TT\ < 0.5), with K = Bn/An- The K values are 
shown in Refs. [5,6]. 

The higher moments of the jet broadenings are ex- 
cluded because of the incomplete [29] description of their 
distributions in general. As in [5] the first moment of i?w 
is excluded because the universality of the fit parameters 
could not be confirmed. 

Thus we combine the results of as(M2o) and 
aoifii) from the moments {{1-T)^), (C^), (B^), 
{iy23V)-{iyf3?), and (A/2), (M^), The parameters are 
consistent within total errors, and the combination proce- 
dure follows that used in [5,6]I3 The correlations between 
as(M2o) and 00(^*7) have rather constant values from - 
0.80 to -0.93 and thus the combination is done individually 
for as(M2o) and ao(/i7). The averages are 

ag(A/z«) = 0.1183 ± 0.0007(stat.) ± 0.0016(exp.) 
±0.0011(had.)+°;°[;^^(^^) 

= 0.1183 ±0.0056(tot.), 
ao(M/) = 0.493 ± 0.006(stat.) ± 0.008(exp.) 

* The covariances of the statistical errors are estimated 
by PYTHIA at hadron level at 91.2 GeV, the covariances 
of the experimental systematic uncertainties are taken as 
Eij = Min{(T^xp.,i , c"oxp.,j}, and hadronisation and renormal- 
isation scale uncertainties are found by repeating the analysis 
with varied parameters. 
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±0.050(had.)_o:oM(2^M) 
= 0.493 ±0.058(tot.). 
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Fig. 3. Measurements of as(Mzo) and ao{^i) from moments 
of six event shape variables at PETRA and LEP energies. The 
inner error bars — where visible — show the statistical errors, the 
outer bars show the total errors. The dashed lines indicate the 
weighted averages, the shaded bands show their errors. Only 
the measurements indicated by solid symbols are used for the 
averages 



4.2 Shape function 

In [8] values for the shape function parameters were de- 
termined by comparing a distribution of Mfj at 91 GeV 
with the NLO-I-NLLA prediction combined with power 
terms from the shape function to be Ai ~ 1.22 GcV, 
Aa = 1.70 GeV2 and 6X2 from ^A2(10GeV) ^ 1.4 (GcV)^ 
to (5A2(100GeV) = 1.2 (GeV)^. 

We perform fits to the moments ((1 — T)^), (C^) and 
(Mg) with two free parameters as{Mzo) and Ai. Fits to 
the higher moments ((1 — T)^), (C^) and (M^) are done 
with an additional parameter. Because of the weak energy 
dependence of SX2 in (j22p wc substitute the numerator 
A2 + 5X2{Q) by the fit parameter A2. The fits to (C^) and 
(Af^) are not sensitive to the parameter A2 while the fit 
to ((1 - T)2) gives A2 =_(1.2 ± 1.3) (GeV)^. For further 
analysis we therefore set A2 = 0. 

Fig. m shows the comparison of the data with the pre- 
diction. The data are fitted well with x^/d.o.f. — based 
on the statistical errors — in the order of one. As before, 
the experimental systematic uncertainties are estimated 
by the minimum overlap assumption, and the renormali- 
sation scale uncertainty by setting = 0.5 and = 2.0 . 




so 100 150 200 SO 100 150 200 

Q[GeV] Q[GeV] 

Fig. 4. Fits of the shape function prediction with A2 = to 
JADE and OPAL measurements of first and second moments 
of 1 — T and C and second and fourth moments of Mh. The 
solid line shows the fitted prediction and the dashed line shows 
the pure NLO contribution. The inner error bars show the sta- 
tistical uncertainties used in the fit and the outer error bars 
show the combined statistical and experimental systematic er- 
rors. Most of the error bars are smaller than the data points 



Tabic [5] and Fig. [5] contain the results for q;s(A/2o) and Ai 
from the standard measurement and the systematic vari- 
ations. 

The values of q;s(A/zo) show the typical increase for 
the higher moments of 1 — T and C. The values of Ai 
are consistent for the first two moments of every event 
shape variable. Furthermore they are consistent for the 
two-hemisphere variables 1—T and C. However, the values 
from Mh are not compatible. As the respective predictions 
are similar, we refer to the high values of q;o(^/) from 
the moments (M^) and (M^) in the dispersive model, see 
Subsect. 14.11 The large power corrections compensate for 
the incomplete perturbative NLO description. Since the 
values are not consistent and no systematic evaluation of 
the theory uncertainty is available, we do not combine any 
of the fit parameters. 
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Fig. 5. Measurements of as(M2o) and Ai from moments of 
three event shape variables at PETRA and LEP energies. The 
inner error bars — where visible — show the statistical errors, the 
outer bars show the total errors 

The Al values from 1 — T, C and, less pronounced, 
from Mh are smaller than the results in [8] . A substantial 
difference of [8] with our work is the inclusion of the NLL A 
approximation in the distribution prediction, which is not 
available for the moments. 



4.3 Variance of event shape variables 

From the first and second moments of event shape vari- 
ables y the variance Var(y) = (j/^) — (y)^ has been cal- 
culated [40]. The experimental uncertainties have been 
determined by systematic variations analogously to the 
moments at hadron level [5,6]. The variance of C, Bt, 
i3w: -^H and becomes larger with rising Q while the 
variance of 1 — T becomes smaller, se^ Fig. [H] and [40]. 
Studying other observables defined in [6] we find that the 
variance of Tmaj. and O becomes larger with rising Q while 
the variance of Tmin., S', -Ml and becomes smaller, 
see [40]. This behaviour is determined by the multi jet re- 
gion of the distribution [41]. The evolution is reproduced 
by Monte Carlo models qualitatively well; however data 
and models differ by typically up to five standard devia- 
tions. 

Fig. [6] shows the comparison of the prediction ([23|) 
using the dispersive model and data of 1 — T, C, B^, 
Bw, ?/23 and M^. The fits of the free parameter as{Mio) 
use central values of the measurements and the statisti- 
cal errors. The energy evolution of prediction and data for 
Var(i3w), Var(l — T), and Var(?/^) coincides qualitatively 
(in case of Var(j/P3) except for the lowest energy points 
of 14 and 22GeV), but not for Var(C), Var(BT) and 
Var(My). The x^/d.o.f. values arc large and the as{Mzo) 
values are not compatible with each other or with estab- 
lished QCD analyses. 

The predictions using the shape function differ only 
for Var(C) from the dispersive model. Fig. [7] shows the 
comparison of this prediction (|24p and data. The fit 
gives xVd.o.f. = 50/8, as(Mzo) = 0.0963 ± 0.0004 and 
Al = (0.824 ± 0.025) GeV (statistical errors). The corre- 
lation of asiM^o) and Ai is 0.80. In the shape function 
description the Q evolution of the perturbative descrip- 
tion deviating from the data is corrected at low energies 
by the term oc Ai . Thus the steep energy evolution at low 
Q is reproduced successfully. However, the fitted values 
as(M2o) and Ai are not compatible with the correspond- 
ing values from (C^) or (C^). 



4.4 Single dressed gluon approximation 

The hadron level data are compared with the five different 
predictions which result from truncating ((28|) after order 
0{a^)...0{a^). The data are described by the predictions 
well. Fig. [8] shows the comparisorF^I of data and predic- 
tion for truncating after a^. The fit results for as and Vi 
arc given in Tabic [51 The coefficients K2...K4 of the more 
strongly suppressed power correction in all instances turn 
out to be compatible with zero. 

For every moment order we study specific properties to 
find the series describing the data best. The x^/d.o.f. val- 
ues from the various expansion orders do not differ signifi- 

^ In case of 1 — T this is not seen very clearly in the plot, 
but significantly in a straight line fit. 

Comparisons for second to sixth expansion order, and more 
discussion, can be found in [41]. 
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Fig. 6. Comparison of the NLO prediction with JADE and OPAL measurements of the variance of 1 — T, C, Bt, -Bw, 2/23 
and M^. The inner error bars show the statistical uncertainties used in the fit and the outer error bars show the combined 
statistical and experimental systematic errors. The line shows the NLO prediction with fitted value of as{M2,o) 



cantly for any moment order and thus do not discriminate 
between different expansion orders. 

((1 — T)^): For the first thrust moment at the order 
"^max = 5 the power correction is compatible with 
zero, then it becomes negativeEB The values of 
as{Mio) and the power correction reach their mini- 
mum at rriyaax = 5, so the convergence appears best at 

((1 — T)^): Including higher perturbative orders the 
asiM^o) and 1^2 values from the second thrust mo- 
ment decrease. The purely perturbative description ap- 
pears continually more accurate, we therefore choose 

((1 - T)3}, ((1 - T)'^): Fitting the third or fourth thrust 
moment, to compensate the negative perturbative 
terms the power correction increases when increasing 
the truncation order. Thus the NLO description ap- 
pears more accurate than the SDG approximation, and 
we choose m^ax = 2. 

Fig. m compares the perturbative terms and the power 
correction terms. In case of the first moment the power 
correction is compatible with zero; the purely perturbative 
prediction from the SDG approximation already describes 
the data well. The power correction is larger for the third 
and fourth moments, but is still small compared to the 
LO and NLO perturbative terms down to Q = 10 GeV. In 
case of the second moment the higher SDG terms up to 
0{aP) are comparable with the power correction. 

There is no distinct sign of the power term. It can be pos- 
itive or negative, or change its sign when using different regu- 
larisations [34]. 



4.5 Measuring as{Mzo) in the SDG approximation 

The as{M-z,o) values steeply rise with moment order in this 
model as observed before in Sect. 14.11 14.21 The problem 
of describing the higher thrust moments consistently — 
i.e. with compatible fit parameters — is not solved in this 
model. We measure the strong coupling in this model only 
from ((l-T)i). 

The central value of the measurement and its statisti- 
cal error are shown in table [S] 5. To de- 
termine the experimental systematic uncertainty, the fit 
is repeated employing the minimum overlap assumption 
for the experimental uncertainties of the data. To esti- 
mate the theoretical uncertainty we study the following 
contributions: 

— In order 0{aP) of the perturbative prediction ((27)) the 
unknown coefficient Pi appears. We set it to the simple 
Fade cstimatcF^ (3^ = /3|//32 instead of (34 = 0. This 
does not change any of the four cited digits of the 
as(Mzo) fit result. 

— Instead of three loops we use the four known loops for 
calculating the Q evolution of as . This changes the fit 
result by less than 0.0001. 

— The 0{a^) coefficient Bi is approximated by SDG 
within 34% [10], and the approximation is expected to 
be complete in high order [9] . Thus we vary the coeffi- 
cient Ci by -1-34% or —34%. This results in changes of 
as(Mzo) by -0.0027 or -^0.0031. This variation is con- 
sistent with truncating after expansion order Wmax = 3 

In this approximation the coefficients (3i follow a simple 
geometric series. 
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Fig. 8. Fits of the SDG prediction including power correction 
to JADE and OPAL measurements of first to fourth moment of 
1 — T. The perturbative part of the prediction has been calcu- 
lated in 0{a^). Superimposed are the perturbative prediction 
in leading, next-to-leading, and maximum order as indicated 
in the figure. The inner error bars show the statistical uncer- 
tainties used in the fit, the outer error bars show the combined 
statistical and experimental systematic errors 



or mmax = 6. A rcnormalisation scale variation is not 
studied and is expected to cause a small effect in a 
prediction of fifth order. 
— As i>i is one of the two fit parameters, the hadronisa- 
tion error is part of the fit error. 

The result is 

as(Mz») = 0.1172 ± 0.0007(stat.) ± 0.0016(exp.) 

+0.0031^-1 N 

= 0.1172 ±0.0036(tot.). 

It is our most precise value. The variation of the NNLO 
coefficient Ci by ±34% is better motivated than the some- 
what arbitrary variation of the rcnormalisation scale fac- 
tor Xfj, = 0.5, 2.0 . It leads to a smaller perturbative un- 
certainty than in the NLO analyses discussed above or 
in [5]. 

In [9] the SDG prediction was fitted to (1 - T) 
data, notably of the PETRA experiments TASSO und 
MARK J at cms energies down to 12 GeV. The results 
were as(Afzo) = 0.110 und Ai =_ 0.62 GeV. In the old 
PETRA measurements e+e^ bb events have not been 
subtracted. They enlarge (1 — T) in an effect decreasing 
with energy faster than \/Q. This results in an increase 
of Al and a decrease of as(A/zo). 



5 Summary 

As an alternative to the description of hadronisation 
by Monte Carlo models we studied analytical models 
of hadronisation. These descriptions require no laborious 
tuning of a Monte Carlo model but fitting of ag and an 
additional parameter. In this way, their consistency can 
be quantified. 

We studied the dispersive model, the shape function 
and the single dressed gluon approximation, including 
the NLO predictions. All models describe the mean val- 
ues of all studied event shape variables at cms ener- 
gies Q = 14...209GeV well. The dispersive model and 
the shape function also describe the higher moments of 
the one-hemisphere variables well. The measurements of 
ag(Mzo) and ao(M/) iii the dispersive model from the 
mean values of 1 — T, C, -Bt, all studied moments of y^^ 
and (M^), (A/^) confirm the universality of these param- 
eters within total errors. 

Higher moments of the two-hemisphere variables 1 — T, 
_Bt and C can only be described by higher values of 
ag (M20 ) , moments of i3w and Mh only by higher values 
of the power correction coefficient. Both deviations appar- 
ently compensate for the deficiencies of the perturbative 
predictions. 

Averaging the parameters from ((1 — T)^), (C^), {B\), 
((2/2°3)')-((2/2°3)')> and {M^), (A/*) in the dispersive 
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Carlo model has only been used for correcting the data 
for bb background, taking into account its uncertainty in 
the experimental systematic error. 

Moments are an interesting alternative to distribu- 
tions: Specific parts of phase space are tested selectively 
and the energy evolution shows up clearly. This evolution 
can be studied thoroughly by the combination of the ex- 
periments JADE and OPAL. As in [5] the most interesting 
energy range is provided by the JADE experiment, as the 
perturbative and non-perturbative effects both scale in- 
versely with the energy. 
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Fig. 9. Q evolution of the employed perturbative terms 
from (|28|) and of the leading power correction in powers of the 
cms energy Q for the case of first to fourth moment of thrust. 
An...F„ are the coefficients of the terms in 1st... 6th order of 
the coupUng a as defined in (|26p . The parameters as(Afzo) and 
Vi are fitted to the data 



model results in 

as(Mzo) = 0.1183 ± 0.0056(tot.) , 
ao(Ai7) = 0.493 ± 0.058(tot.) . 

Event shape moments are barely sensitive to more strongly 
suppressed power corrections resulting from the shape 
function. In general, a detailed understanding of non- 
perturbative effects cannot reach any further than the reli- 
ability of the underlying perturbative structure which was 
studied in [5, 6]. 

For some variables the variance shows an unexpected 
energy evolution. Of the variables analysed here only the 
1 — T distributions are broader at low cms energy; the 
M^, C, iJx and i?vv distributions are narrower at low 
energy. The energy evolution of the variances is not given 
correctly by the pure NLO prediction. The description by 
an additional power correction from the shape function is 
not quantitatively correct. 

The single dressed gluon approximation shows the 
same problems describing the higher thrust moments as 
the other models. Fitting the prediction for ((1 — T)^) in- 
cluding terms up to fifth order gives power corrections 
compatible with zero and 

as(Mz") = 0.1172 ± 0.0036(tot.) . 

With a competitive total error of 3% this is our most pre- 
cise as measurement. It agrees with the world average of 
as(Mzo) = 0.1189 ± 0.0010 [42]. For this result a Monte 
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Table 3. Measurements of Qs(A'f2o) and Qo(/ii) in the dispersive model from first moments of six event shape variables over 
the full analysed range of cm. energy, 14...209GeV. The xVd-o.f. values are based on the statistical errors only, see text for 
further discussion 





((i-r)i) 










as(Mzo) 


0.1234 


0.1228 


0.1200 


0.1142 


0.1181 


Statistical error 


0.0006 


0.0004 


0.0005 


0.0006 


0.0004 


Experimental syst. 


0.0015 


0.0012 


0.0011 


0.0015 


0.0016 


variation: = 2.0 
— 0.5 


+0.0073 
-0.0062 


+0.0073 
-0.0063 


+0.0047 
-0.0041 


+0.0004 
+0.0009 


+0.0050 
-0.0038 


fii variation; = 1 GeV 
fii =3 GeV 


+0.0025 
-0.0019 


+0.0030 
-0.0022 


+0.0017 
-0.0013 


+0.0008 
-0.0007 




M variation: M - 20% 
M + 20% 


+0.0011 
-0.0011 


+0.0014 
-0.0013 


+0.0008 
-0.0007 


+0.0004 
-0.0004 


— 


Theoretical syst. 


0.0078 


0.0080 


0.0051 


0.0012 


0.0050 


aoifJ-i) 


0.490 


0.433 


0.514 


0.604 




Statistical error 


0.007 


0.003 


0.009 


0.015 




Experimental syst. 


0.010 


0.005 


0.009 


0.014 




Xfj, variation: x^ = 2.0 
Xfj, = 0.5 


-0.022 
+0.031 


-0.024 
+0.033 


+0.013 
+0.011 


-0.049 
+0.147 




M variation: M - 20% 
M + 20% 


+0.047 
-0.033 


+0.034 
-0.025 


+0.053 
-0.037 


+0.086 
-0.057 




Theoretical syst. 


0.056 


0.048 


0.055 


0.170 




Correlation asiM^^o), ao{^ii) 


-0.86 


-0.80 


-0.93 


-0.90 




xVd.o.i. 


47.8/10 


66.3/10 


106/10 


75.3/10 


37.7/10 
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Table 4. Measurements of qs(M2o) and ao{i-ii) in the dispersive model from second to fifth moments of six event shape 
variables over the full analysed range of cm. energy, 14...209GeV. The xVd-o.f. values are based on the statistical errors only, 
see text for further discussion 













(MA) 




as(Mzo) 


0.1182 


0.1158 


0.1173 


0.1160 


0.1154 


0.1144 


Statistical error 


0.0009 


0.0007 


0.0013 


0.0017 


0.0010 


0.0020 


Experimental syst. 


0.0021 


0.0024 


0.0020 


0.0025 


0.0034 


0.0033 


variation: — 2.0 
— 0.5 


+0.0054 
-0.0042 


+0.0042 
-0.0033 


+0.0053 
-0.0041 


+0.0050 
-0.0038 


+0.0051 
-0.0042 


+0.0047 
-0.0036 


fii variation: /i/ = 1 GeV 
^7=3 GeV 




+0.0013 
-0.0010 






+0.0011 
-0.0009 




M variation: M - 20% 
M + 20% 


— 


+0.0006 
-0.0006 


— 


— 


+0.0005 
-0.0005 


— 


Theoretical syst. 


0.0054 


0.0045 


0.0053 


0.0050 


0.0053 


0.0047 


ao(fJ-i) 




0.601 






0.579 




Statistical error 




0.011 






0.015 




Experimental syst. 




0.012 






0.022 




Xfj, variation: a;^ — 2.0 
Xfj, = 0.5 




-0.022 
+0.037 






-0.022 
+0.034 




M variation: M - 20% 
M + 20% 




+0.083 
-0.056 






+0.079 
-0.054 




Theoretical syst. 




0.091 






0.086 




Correlation as(Afzo), 00(^7) 




-0.88 






-0.82 




xVd.o.f. 


18.3/10 


25.2/10 


21.0/10 


18.7/10 


26.4/10 


15.2/10 
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Table 5. Measurements of ots{M,^a) and Ai in the shape func- 
tion model from the first and second moments of three event 
shape variables (second and fourth moment in case of Mh) over 
the full analysed range of cm. energy, 14. ..209 GeV 







((1 - TV) 


\^ / 




as(Mzo) 




0.1304 


0.1305 


0.1193 


Statistical error 




0.0008 


0.0006 


0.0008 


Experimental syst. 




0.0018 


0.0015 


0.0027 


variation: = 


2.0 
0.5 


+0.0069 
-0.0054 


+0.0064 
-0.0050 


+0.0038 
-0.0025 


Ai[GeV] 




0.499 


0.393 


1.010 


Statistical error 




0.034 


0.024 


0.047 


Experimental syst. 




0.065 


0.045 


0.075 


Xfj, variation: = 

X ^ 


2.0 
0.5 


+0.053 
-0.067 


+0.063 
-0.075 


+0.044 
-0.045 


Correlation asiM^si), Ai 


-0.96 


-0.96 


-0.94 


xVd.o.f. 




23.3/8 


22.6/8 


19.6/8 








1 




as(Mzo) 




0.1426 


0.1456 


0.1184 


Statistical error 




0.0014 


0.0008 


0.0012 


Experimental syst. 




0.0028 


0.0023 


0.0041 


x^ variation: = 
Xfi — 


2.0 
0.5 


+0.0111 
-0.0091 


+0.0117 
-0.0097 


+0.0049 
-0.0037 


Ai[GeV] 




0.516 


0.364 


0.964 


Statistical error 




0.056 


0.023 


0.079 


Experimental syst. 




0.096 


0.052 


0.133 


Xf_i variation: = 

X f^i — 


2.0 
0.5 


+0.044 
-0.052 


+0.063 
-0.054 


+0.031 
-0.030 


Correlation as{M2 


;o), -^1 


-0.95 


-0.91 


-0.92 


xVd.o.f. 




20.9/8 


18.3/8 


20.2/8 
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Table 6. Measurements of Qs(Af2o) and Un from first to fourth moment of 1 — T over the full analysed range of cms energy, 
14. ..209 GeV. The single dressed gluon approximation is used in maximum order mmax. The errors are statistical. Units of Un 
are GeV torn = 1 and (GeV)^ for n = 2, 3, 4 













correlation 


/d o f 










as(Mzo), Vn 


((1 


-Tf) 


2 


0.1314±0.0008 


0.430±0.036 


-0.96 


29.1/10 






3 


0.1209±0.0007 


0.352±0.039 


-0.97 


30.0/10 






4 


0.1178±0.0007 


0.223±0.045 


-0.97 


29.3/10 






5 


0.1172±0.0007 


0.035±0.058 


-0.98 


26.4/10 






6 


0.1182±0.0001 


-0.293±0.010 


0.88 


24.0/10 


((1 


~Tf) 


2 


0.1463±0.0009 


1.855±0.455 


0.73 


23.8/10 






3 


0.1426±0.0009 


1.762±0.454 


0.74 


24.0/10 






4 


0.1416±0.0008 


1.644±0.458 


0.73 


24.2/10 






5 


0.1412±0.0008 


1.546±0.460 


0.73 


24.3/10 






6 


0.1410±0.0008 


1.465±0.470 


0.73 


24.4/10 


((1 


-Tf) 


2 


0.1523±0.0015 


0.273±0.136 


0.73 


20.7/10 






3 


0.1525±0.0015 


0.274±0.137 


0.73 


20.7/10 






4 


0.1531±0.0015 


0.286±0.137 


0.73 


20.7/10 






5 


0.1533±0.0015 


0.298±0.136 


0.73 


20.6/10 






6 


0.1534±0.0015 


0.307±0.135 


0.74 


20.6/10 


((1 




2 


0.1581±0.0017 


0.059±0.024 


0.55 


17.6/10 






3 


0.1625±0.0018 


0.063±0.024 


0.55 


17.5/10 






4 


0.1645±0.0019 


0.075±0.023 


0.56 


17.4/10 






5 


0.1652±0.0019 


0.085±0.022 


0.58 


17.3/10 






6 


0.1655±0.0019 


0.095±0.021 


0.60 


17.3/10 
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